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Abstract 

UVES and HIRES high-resolution spectra of comet 9P/Tempel 1 are used to investigate the 
impact and rotational light curves of various species with a view toward building a simple model 
of the distribution and activity of the sources. The emission by OH, NH, CN, C3, CH, C2, 
NH2 and 01, are analyzed, as well as the light scattered by the dust. It is found that a simple 
model reproduces fairly well the impact lightcurves of all species combining the production of 
the observed molecules and the expansion of the material throughout the slit. The impact 
light curves are consistent with velocities of 400-600 m/s. Their modelling requires a three-step 
dissociation sequence "Grand-Parent — > Parent — > Daughter" to produce the observed molecules. 
The rotational light curve for each species is explained in terms of a single model with three 
sources. The dust component can however not easily be explained that way. 

Keywords: Deep Impact — Comets — 9P/Tempel 1 
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1 Introduction 

The Deep Impact (DI) experiment was a first in cometary astrophysics. Instead of observing 
objects in a relatively steady state, it was possible to watch the instantaneous release of a 
large quantity of cometary material and its subsequent evolution in timescales of hours and 
days (A'Hearn et al. 2005). The fact that the time of origin of the observed effects is known, 
provided an extraordininary advantage in studying the complex evolutive mechanisms governing 
the physics and chemistry in the coma. The international collaboration also provided the unique 
opportunity to monitor a comet during weeks with an armada of telescopes (Meech et al. 2005). 
Here, we present the analysis of data obtained over 15 nights spead before, during and after 
the Deep Impact event (2005 July 4, 05:52 UT as seen from Earth) with the high-resolution 
spectrographs of the ESO VLT (European Southern Observatory, Very Large Telescope, Cerro 
Paranal, Chile) and Keck (Hawaii) observatories. We investigate the impact and the rotational 
light curves of various emissions and provide a simple model for the activity of the sources. The 
emission by OH, NH, CN, C3, CH, C2, NH2 and OI, are analyzed, as well as the light scattered 
by the dust. 

The observations are presented in Section |2J The light curves observed immediately after the 
impact at the Keck I telescope are presented in Section |3] together with their interpretation by 
a simple model. In Section |3J these light curves are again used to build a model also for the 
rotational light curves that resulted from the data collected during 13 nights at the VLT. Finally 
we summarize our conclusions in Section [5J 
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2 Observations and data reduction 

The high-resolution spectra of comet 9P/Tempel 1 obtained with the VLT UVES (ESO UltraVi- 
olet Echelle Spectrograph) and Keck I HIRES (High Resolution Echelle Spectrograph) between 
May and July 2005 are described in Jehin et al (2006), and in Tables ^ and EJ The slits used 
at both telescopes had similar areas (10" by 0.44" for UVES, 7" by 0.86" for HIRES) and were 
generally centered on the nucleus, which allows some comparison of the photometric light curves. 
The reductions were made in the usual manner for cometary spectra. The echelle package of 
IRAF (Image Reduction and Analysis Facility) and the UVES pipeline (vl.4) implemented in 
the ESO-MIDAS (Munich Image Data Analysis System) software package (Ballester et al. 2000) 
were used to extract the various orders. The dust and twilight contributions were subtracted 
after proper scaling and Doppler shifting of an appropriate solar (UVES) or solar analog (HIRES) 
spectrum. This was done order by order so that linear wavelength shifts and simple scaling 
factors could be used. The latter was derived in an iterative procedure until nulling of all 
remaining solar features. 

The spectrographs were not designed to be used as high precision photometers, and the observa- 
tions were not conducted with the photometric issue as a priority. The centering on the nucleus 
was not always perfect and the slit width was only a fraction of an arc second. Some UVES 
spectra were intentionally offset from the nucleus by between 1 and 2 arcsec (see Table This 
can lead to significant effects on components with a sharp spatial profile (dust). The position of 
the nucleus along the slit was not quite stable in the HIRES data. A slight drift along the slit 
occured during the impact night, partly because of the differential atmospheric refraction (the 
slit was aligned along the parallactic angle and there was no refraction corrector). At high air 
mass, in the bluest orders of the last spectra, the nucleus was almost at the edge of the slit. 
The main problem was thus to choose the proper extraction aperture. Calculations made with 
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different apertures showed differences of up to 15 % in the photometric light curve during the 
impact night. After some testing we found that the most sensible solution was to use the full slit 
length on both spectrographs, keeping in mind the intrinsic limitations of such measurements. 
This full aperture could be easily defined from the flat frames and it yielded the maximum 
signal-to-noise ratio, which is certainly a determining factor for such a faint comet. 
From these spectra we derived the flux of major emission bands and of the dust scattered solar 
spectrum, assuming a clear atmosphere, and stable detectors. Eight gaseous species were mea- 
sured and analyzed : OH (0-0) at 308 nm, NH (0-0) at 336 nm, CN (0-0) at 388 nm, C 3 complex 
at 405 and 407 nm, CH (0-0) at 430 nm, C 2 (0-0) band head at 516 nm, various lines of NH 2 
at 519, 531, 569, 574 nm, and the green and red oxygen forbidden lines ([OI]) at 558 nm and 
630 nm (UVES only). Table details the wavelength ranges of interest where major emission 
features and the corresponding comparison zones (not always emission free) can be found. The 
wavelengths are given in the rest frame of the nucleus. For OH, NH and CN the integration 
was done on a series of individual lines. The [O I] lines have been deblended from the telluric 
features. Plots of a typical spectrum are shown on Fig. ^ 

UVES data were obtained over 13 nights (10 consecutive nights encircling the impact date). 
This provided complete coverage of the comet's light curve over several rotational periods. As 
the CN band was covered in all UVES setups, it is the species with the most observations. Other 
molecules have less dense coverage in the phase diagram. After correcting for the higher efficiency 
of new mirror coatings and for the strong Swings effect closer to the comet's perihelion, we could 
directly compare CN fluxes measured in June and July 2005. The comet activity appeared to 
be slightly higher (about 3 %) in June, for which we corrected the CN data. However, for the 
other species we used only the July observations, as we could not do reliable adjustments owing 
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to less available data and/or the lack of accurate models. 

The dust contribution has been accurately evaluated in each spectrum in order to perform 
an accurate continuum subtraction. However, as the small slit width of UVES lead to large 
variations when the slit was not perfectly centered on the nucleus and/or when the seeing 
changed, we kept only the well-centered spectra (see Table ^) for evaluation. The HIRES data 
were obtained with a wider slit and therefore did not show this effect. The main HIRES data 
set covers the four hours following the impact (see Table 3). In addition three HIRES spectra 
obtained on May 30 are available . 

The relative variations in the brightness of the dust appear to be wavelength independent over 
the whole spectral range, both in the UVES and HIRES impact data. The data presented here 
are the same as in Jehin et al (2006). 

3 Modelling the impact light curve 

The impact light curves were deduced from the HIRES and UVES spectra obtained during the 
hours (or nights) following the impact. They had to be corrected for the natural rotational 
variations, which were most important for the spectra obtained during the first hours after the 
impact as they coincide with the very phase when the natural light curves exhibit the fastest 
variations of the 1.7 d cycle. The resulting impact light curves are illustrated in Figs. 121 4101 

3.1 A simple model 

We assume that the observed species -the daughter- originates from a single parent which itself 

comes from a single grandparent in a continouous sequence of dissociation (the three species 

coexist at any t). The grandparent appears necessary to introduce a delay in the light curves 
1 Those data are publicly available at http://msc.caltech.edu/deepimpact/ 
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(cf. Sect. 3.2). Within this framework, we may write 

dn g = —agTigdt 

drip = +a gp n g dt — a p n p dt (1) 
dn d = +a pd n p dt — a d n d dt 

from which we derive 

n d a pd a gp e~ adt H apd Qgp e~ apt H apd a§p e - " 9 * (2) 

n go (a p - a d ){a g - a d ) (a d - a p )(a g - a p ) (a d - a g )(a p - a g ) 

In these equations, n is the number of molecules, the subscripts d, p, g standing for daughter, 
parent and grandparent respectively. n go is the number of grandparents produced at once by 
the impact. The grandparent can be anything, e.g. icy grains, complex or simple molecules, as 
long as it obeys Eq. 1. a = 1/r are inverse lifetimes; a g , a p and a d characterize the destruction 
or dissociation of the grandparent, parent and daughter; a gp and a pd refer to the production 
of the parent from the dissociation of the grandparent and to the production of the daughter 
from the dissociation of the parent, accounting for the possibility that the grandparent (resp. 
parent) may dissociate into several parents (resp. daughters). We then have a g = J2 P a gp an d 
a p = J2d a pd- If the dissociation lifetime of the daughter is very long compared to the timescale 
of the event (a d ~ 0) and if the grandparent does not exist or has a very small lifetime, Eq. 2 
will simplifly to 

= 1 - e -«pt (3) 

which is useful to get a first insight into the observed light curves. In this the number 

of parents produced by the impact. 

Eqs. 2 and 3 are only valid if all the expanding material is within the slit area. When a part 
of the material moving at velocity v exits the slit width w, i.e., at t > t w = w/2v, the number 
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of particles within the slit decreases with time. Assuming that the expanding material forms 
a shell of radius v t, the volume inside and outside the slit area may be readily computed such 
that for t > t w 

n d = n d— (4) 

where is given by Eq. 2 or 3 (when t < t w , n' d = n^). Although established for a spherical 

expansion, this relation is also valid for instance for a hemispherical flow. When t > ti, the 

material fills both the slit width w and the slit length /. Assuming / » w, we may write 

' t w 2 ti 

n d = n d — -arcsm(-) (5) 

t TT t 

which becomes 

' 2 t w ti 

n d = n d - -p (6) 

when t S> t\. In the case of the HIRES data, ti = 8.1 t w . Assuming the medium to be optically 
thin, the observed flux is simply proportional to n' d (t) integrated over the exposure time. 
Since it is unlikely that all the material is moving at a single velocity v, we introduce a velocity 
distribution. For convenience we use 

dn go = NgQ -^e-^ 2 ^dv (7) 

where the mean velocity v and its standard deviation a v are parameters which can be estimated 
from the observations. It is important to note that in this simplified model we assume that the 
daughter, parent and grandparent species expand with the same velocity. 
To summarize, we have 

n d (t) = j N go BD(v)R(t)S(v,t)dv (8) 
integrated over the velocity distribution, with 

D{v) = —^-= e-^fl-l ; ( 9 ) 
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R(t) 



(a p - a d )(a g - a d ) 



e -a d t + 



C^p Gt g 



(a d - a p )(a 6 



_ e -<* P t + 



dp OLg 



-a a t 



(a d - a g )(a p - a g ) 



(10) 



and 



S(v, t) = 1 when t w > t 

S(v,t) = — when t w < t and ti > t (11) 

S(v,t) = — — arcsin( — ) when t[ <t 
t tt t 

where t w = w/2v, ti = l/2v and B = a pd a gp / 'a p a g . Free parameters are v, a v , a d , a p , a g and 
N go B. Since absolute fluxes are not measured, N go B is just a scaling factor. As seen in Eq. 10, 
R(t) is identical if we exchange the values of a p and a g , such that parameter fitting can be done 
with either a g /a p > 1 or a 9 /a p < 1. 

3.2 Interpretation of the light curves 

CN — Fig. |2 illustrates the observed post-impact light curve of CN (hereafter the daughter 
molecule). It is characterized by a rather steep rise and a smooth decline. The maximum 
corresponds to the time at which the material fills either the slit width or the slit length. At 
times t > t w and t > tj, the light curve can be seen as a combination of the production of 
CN (Eq. 2) and the dilution of the expanding material beyond the slit area (Eqs. 4 or 5). At 
large t we may neglect, to first approximation, the effect of the grandparent as well as the CN 
dissociation lifetime, which is very long (Fray et al. 2005 and references therein) — probably 
longer than 4 days (Keller et al. 2005). Therefore Eq. 3 holds at large t. In this case, the 
number of daughter molecules may be constant after some time (as long as a p t <C 1 in Eq. 3) 
or proportional to t (if a p t 1). We then expect the shape of the decline of the light curve 
to be roughly intermediate between a constant and i" 1 if the bump corresponds to the filling of 
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the slit width (Eqs. 3 & 4) or between t -1 and t~ 2 if the bump corresponds to the filling of the 
slit length (Eqs. 3 & 6). However, since the velocity of the flow is not unique the situation is 
less cut-clear and an independent knowledge of the bulk velocity of the flow is needed. 
In the case of CN, the analysis of the spatial profiles has led us (Jehin et al. 2006) to estimate the 
velocity of the expanding material to v — 400 m/s (with a range of projected velocites from 250 
to 650 m/s). Since the 0.86" x 7.0" HIRES slit corresponds to 562 km x 4570 km at the comet 
at the epoch of the observations, the slit length is filled after ti = 0.066 d, which corresponds 
to the bump seen in the light curve. The fact that the decline of the light curve is not as steep 
as t~ 2 indicates that the CN production still continues at t > t\ which requires a rather long 
lifetime for its parent. If the material is expanding at a single velocity (a v we will expect 

a sharp decline just after filling the slit, i.e., a relatively narrow bump. The fact that such a 
bump is not observed requires some smoothing, i.e. a large velocity dispersion. After some trial 
we find that a v > 100 m/s is needed to reproduce the data, which is in agreement with initial 
estimates (Jehin et al. 2006). 

As illustrated in Fig. [21 modelling the CN light curve with only a parent molecule results in 
a much faster rise after the impact and a clear discontinuity at t ~ t w , when the material 
leaves the slit width. Such models cannot reproduce the data whatever values are adopted for 
the lifetimes. A grandparent must be added to introduce a delay in the light curve just after 
the impact. In Fig. [2^, we show reasonable fits to the data with different combinations of a p 
and a g /a p , indicating some degeneracy between a p and a g . The value of the parent lifetime 
needed to reproduce the data, does not only depend on the ratio a g /a p , but also on the adopted 
velocity distribution. Fig. shows an equally acceptable fit with a v = 200 m/s. In general, 
adding material at lower velocities enhances the number of molecules in the slit at t > ti, and 
consequently requires a shorter parent lifetime to fit the data. 
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Although we are able to reproduce fairly well the observed light curve, the parent lifetime 
t p = l/a p that we derive appears very uncertain owing to the strong dependence on both 
the velocity distribution and the ratio a g /a p . Without additional data and/or independent 
constraints, the CN parent dissociation lifetimes, whcih we estimate from acceptable fits, are 
roughly in the range 0.1 - 10 d. Fortunately, additional measurements from UVES observa- 
tions were performed several hours after the impact allowing us to constrain further the model 
(Fig. 12b). We find that models with 0.2 d" 1 < Q p < 4 d 1 are favoured, which corresponds to 
CN parent dissociation lifetimes of 2 10 4 s < r p < 5 10 5 s. Even with this additional constraint, 
the derived lifetimes are accurate only within orders of magnitude. They are consistent with 
other estimates (Keller et al. 2005, Fray et al. 2005) and with the hypothesis that HCN is the 
parent of CN. The CN grandparent destruction lifetime is even more uncertain, providing no 
information on its nature. 

Figures 

EH2l illustrate the light curves of the other molecules and [OI]. The modelling is essen- 
tially similar to the one discussed in detail for CN. When necessary, the velocity distribution is 
fine tuned in order to reach the best fit. Indeed, differences can arise in, e.g., the dissociation 
process. Additional details are briefly described below while the derived parameters are gathered 
in Table H 

OH — In Fig. |31 the photometry obtained with the Optical Monitor (OM) of the XMM-Newton 
satellite (Schulz et al. 2006) has been included. The simultaneous fit of HIRES and XMM-OM 
data requires a velocity of the OH radicals comparable to the velocity of CN. We then adopt 400 
m/s for the material to exit the HIRES slit length at ti = 0.066 d and the XMM-OM aperture 
at t w = 0.056 d. Models of the light curve are not very sensitive to the long OH dissociation 
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lifetime (jd ~ 2 10 5 s at = 1.5 AU, from Morgenthaler et al. 2001) and require both a grand- 
parent and a parent (Fig. I3J . Acceptable models reproduce fairly well the data with the parent 
and grandparent destruction lifetimes given in Table |U which, as in the case of CN, appear to 
be poorly determined. The OH parent dissociation lifetimes we derive largely encompass the 
dissociation lifetime of H2O (~ 8 10 4 s), a possible parent of OH. 

NH2 — The relatively narrow bump in the NH2 light curve (Fig. @J suggests a slightly higher 
velocity with a not too high dispersion. We adopt v = 500 m/s and a v = 100 m/s. A range of 
acceptable combinations of a p and a g is given in Table |IJ The long NH2 dissociation lifetime 
has little effect on the models (r^ > 1.5 10 5 s at = 1.5 AU, from Cochran et al. 1992). The 
derived NH2 parent lifetimes are poorly determined and are consistent with other estimates (e.g. 
Wyckoff et al. 1988, Fink et al. 1991). 

NH — A slightly higher velocity is needed to reproduce the data; we adopt v = 550 m/s and 
a v = 150 m/s (Fig. EJ). The long NH dissociation lifetime has little effect on the models (r^ > 1 
10 5 s from Schleicher et al. 1989). The derived NH parent dissociation lifetimes are consistent 
with other estimates (e.g. Schleicher et al. 1989). 

C3 — Keeping in mind the large uncertainties of the flux measurements (cf. the measurements 
just after the impact), acceptable models can reproduce the data with v ~ 400 m/s (Fig. |SJ). 
A range of acceptable combinations of a p and a g is given in Tabled The long C3 dissociation 
lifetime has little affect on the models (r^ > 3 10 5 s at = 1.5 AU from Cochran et al. 1992). 
The rather short C3 parent dissociation lifetimes we derive are consistent with other estimates 
(e.g., Rauer et al. 2003, Randall et al. 1992). 
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C2 — A higher velocity is needed to reproduce the data; we adopt v = 650 m/s and a v = 200 m/s 
(Fig. Q). However the measurements are very uncertain and should be seen with caution. The 
long C2 dissociation lifetime has little effect on the models (r^ > 3 10 5 s at = 1.5 AU from 
Cochran et al. 1992). The derived C2 parent dissociation lifetimes are consistent with other 
estimates (Combi and Fink 1997 and references therein). 

CH — While the range of the adopted CH parent and grandparent lifetimes (Tabled} reproduces 
reasonably well the major part of the light curve (given the uncertainties) and agrees with other 
estimates (e.g. Cochran et al. 1992), the relatively short CH dissociation lifetime used in the 
modelling (r^ > 1 10 4 s at = 1.5 AU, from Cochran et al. 1992) might be underestimated 
since the residual flux at i — 0.85 d cannot be fitted by the model (Fig. [SJa). However the 
measurements for CH are very uncertain and should be regarded with caution. 

OI — Photodissociation of molecules like H2O, OH, CO2 or CO can produce 01 directly in 
the excited X D and X S states (Festou and Feldman 1981, Morgenthaler et al. 2001). These 
metastable states have a lifetime of ~ 110 s and ~ 0.74 s, respectively, much shorter than the 
photodissociation lifetime of the parent molecules. Those atoms travel only a short distance 
before decaying to the ground state and emitting photons at 6300A (and 6364 A) and at 5577A 
respectively. 

In the modelling procedure, we first considered OI atoms as the daughter species. However, 
given the very short lifetime of the metastable levels, the light curve actually traces the parent 
molecule, such that one stage in the 3-step dissociation sequence used up to now is lost. Such 
models provide poor fits essentially because the computed light curves are rising too fast just 
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after the impact. When considering either OH or CO as the "daughter" species -with rather 
long dissocation lifetimes- and H2O or CO2 as the "parents", coming themselves from unknow 
grandparents released by the impact, acceptable fits of the light curve are obtained (Fig. [HI) with 
various combinations of a p and a g (Table0J). While the information the derived lifetimes provide 
on the parent species is only limited, the fact that a 3-step dissociation sequence is needed to 
reproduce the data suggests that 01 ( 1 S) atoms are not produced directly from the dissociation 
of H2O or CO2, but through OH or CO. Moreover, the 01 light curve is not similar to the one 
of OH (Fig. |HJ). This might indicate a significant contribution from the CO2 — > CO dissociation 
sequence. 

3.3 The dust light curve 

The variations of the dust after the impact are much faster than for the gas, with a very steep 
brightening fFig. I1U|). The maximum is already reached at about UT 06:18 db lOmin, the dust 
being at that time enhanced by a factor of ~ 8.5 with respect to the pre-impact spectrum. The 
decline of the dust emission may also be interpreted as the escape of dust from the slit. This 
would give a projected dust velocity of ~ 0.18± 0.05 km/s, as, in this case, the escape from the 
slit width is the dominant factor. However the slow and quasi linear decline would require a 
rather broad range of velocities - slower than the gas component - and/or complex processes like 
destruction of highly reflective icy grains by sunlight. Similar dust outflow projected velocities 
were deduced from the analysis of the radial profiles (~ 0.13± 0.03 km/s, Jehin et al. 2006) 
and from the evolution of the dust cloud (~ 0.18 — 0.2 km/s, Meech et al. 2005 and references 
therein; ~ 0.13 - 0.23 km/s Schleicher et al. 2006), . 
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4 Modelling the rotational light curve 

A preliminary analysis of the CN light curve from the UVES data (Fig. lllj) revealed a periodicity 
of 1.709 ± 0.009 d and the presence of three main sources or active regions (Jehin et al. 2006). 
We constructed phase diagrams with the above period using the time origin of the impact as 
seen from Earth, July 4, UT 05:52. In order to allow comparisons, all data were normalized to 
unity at impact time, i.e., at phase 0. The spectrum obtained at impact is strongly affected by 
twilight and yields noisier data. Hence, the May data were included in the phase diagram to fix 
the zero point. 

We want to test the hypothesis that the light curve f s (t) of the comet in the wavelength of 
different species (s) can be represented as the response of the comet to a continuous succession of 
infinitesimal impact-like events. This is merely a way of describing the instantaneous production 
rates of the active regions under the influence of the solar radiation modulated by the rotation of 
the nucleus. We assume a linear behavior, i.e., f s (t) is the linear superposition of the transitory 
responses to the series of micro-events: 

/,(*) = /* g s {t')n s (t-t')dt' (12) 

where n s (t) is the impact light curve (Eq. 8) corresponding to species s and g s (t) represents the 
total instantaneous activity of all sources. g s (t) = g s (t + P) and f s (t) = f s (t + P) are periodic 
functions of period P. We assume that there is a finite number N of sources. Their contribution 
can be described as 

N 

ft(0=EM*' -t*,i)+9;0 (13) 

i=l 

where the b S) i(t) = b s ^{t + P) are the rates of activity of the N sources, and g S) o is a possible 
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background of activity. 

The decomposition into discrete sources arises from physical considerations, i.e., the idea that 
the cometary material comes from several areas or vents. It also allows us to use simple functions 
b(t) since one might expect the activity of each source to exhibit periodic rise and fall cycles. 
One may also expect the various sources to have similar behavior, so the same function could 
be used with some phase shift. 

It must be emphasized that, because impact and rotational light curves overlap for several days, 
they must be disentangled through an iterative process, which adds some complexity to the 
modelling procedure. 

The following conclusions were reached using the simplest satisfactory model, i.e., the one with 
fewest free parameters. 

• The two main maxima of the light curves suggest the presence of at least two sources. 
In fact, three sources are necessary to explain the general shape of the light curves, i.e., 
N = 3. Using only two sources leads to a pronounced gap in the light curve. 

• There is no residual activity g s $ = 

• The activity of each source follows the same pattern during the rotation of the nucleus, 
i.e., we can write : 

b s ,i(t) = asM*) ( 14 ) 

with the b s (t) independent of i. 

• This activity of each source can be approximated by an intermittent function, a variant of 
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a triangular wave. With < t Sj i < r Sj 2 < P, we write in the interval < t < P : 

= V r s,l if < t < t Sj1 

b,(t) = { = ( Ts 2 _ t)/( Ts , 2 - T 8 ,i) if r S)1 < t < r S)2 (15) 
= if r Sj 2 < t < P 

This function repeats itself at each interval kP <£<(& + 1)-P- 

• The onset times of each source are the same for each species, t Sj i = U = (0.18P, 0.47P, 0.85^)+ 
kP 

• The relative intensity of each source is the same for all species, i.e., a s j = c s di with 
d{ = (0.4,0.7,1.0). This is equivalent with the relative production rates of each species 
being the same for each source, i.e., the material is the same. 

• The long-term characteristic time-scale of the decay of n s (t) is strongly correlated with 
the overall amplitude of the light curve s. This provides an independent way of predicting 
the residuals in the post-impact nights. 

• Finally, r Sj i and r Sj 2, describing the activity of the sources, are very similar for all species. 
We found that the best fits are obtained with two sets of parameters. For CN, CH and 
[OI] we find r,,i = 0.21P and r Sj2 = 0.39P while for OH, NH, NH 2 , C 2 , C 3 , t sA = 0.19P 
and Tg.2 = 0.35P. The rise and fall are quasi symmetrical. 

The above conclusions apply to the gas only. The dust cannot be described by the same model. 
The phase diagram shows almost no hint of periodicity whereas a fit with the above parameters 
and the response ?\iust(*) observed after the impact would give very strong variations (see 
Fig. 111! bottom). This means that either (i) the temporal variation of the dust luminosity 
observed after the impact, ra^ ust (t), is not representative of the response of the dust ejected 
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from the sources, or (ii) the modulation is masked by a very large residual component Q/dugt o) 
of non- volatile dust. 

5 Discussion and conclusions 

The high-resolution spectra obtained with the UVES and HIRES spectrographs have provided 
photometric data of unexpectedly good quality, in view of the faintness of the target, the limited 
area of the slits and the less than perfect centering. This allowed us to model the impact as well 
as the rotational light curve of this comet which led to the following conclusions. 

• Within the uncertainties a simple model reproduces fairly well the impact light curves of 
all species. The light curve may reflect the combination of the production of the observed 
species and the expansion of the material throughout the slit. 

• All impact light curves are consistent with v ~ 400-600 m/s (a v = 100-200 m/s). NH and 
NH2 have similar light curves and slightly higher velocities, possibly indicating a common 
origin. 

• With these expansion velocities, a three-step dissociation sequence "Grand-Parent — > Par- 
ent — ► Daughter" is necessary and sufficient to model all light curves. 

• The exact nature of the grandparents is unknown. It can be anything which desintegrates 
(or sublimates) as e~ t l 1 ' 9 . The order of magnitude of its destruction lifetime is typically 
10 4 s. It could be common to all species (e.g., icy dust grains released by the impact). 

• Parent lifetimes can only be determined to better than an order of magnitude. They 
largely encompass the various parent lifetimes proposed in the literature. 
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• The gas activity is generated by three sources, whose duration is less than about 0.4 of 
the period, i.e., compatible with a daylight excitation. 

• The activity of each source follows a similar pattern. 

• The progressive increase and decrease of the activity of each source is compatible with the 
modulation of the incidence angle of the solar radiation. 

• The dust observed during the impact disappears rapidly and is probably mostly sublimat- 
ing water ice in fragments with a diverse size distribution. 

• The dust observed after the impact is dominated by another, less volatile component, which 
may consist of the core of slow-moving dust particles that have lost their icy component 
and are weakly bound to the nucleus. 

• The material releasing CN, CH and 01 seems to be present during a slightly longer period 
than that responsible of OH, NH, NH 2 , C 2 , C 3 . 
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Table 1: UVES data observational circumstances 
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r is the heliocentric distance in astronomical units (AU) and A is the geocentric distance, r is the 
heliocentric radial velocity, Exp. is the the total exposure time in seconds, Airm is the air-mass at mid 
exposure, Slit PA is the slit position angle (North to East), Off. is the offset of the slit center from the 
nucleus, PsAng is the position angle of the Sun-Comet vector. The Setup column defines the central 
wavelength of the UVES blue and red arm. 
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Table 2: HIRES data observational circumstances 
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Table 3: The gaseous species analyzed 



Species Central windows Side windows Extinction 
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Table 4: Results of model fitting 
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Figure 1: Spectral features used to measure the various species. They are extracted from 

spectra taken with UVES the night following the impact. The flux is given in arbitrary units. 

The cometary [01] is the small bump on the red wing of the strong telluric line. Most of the 
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band blueward of 517 nm is C2, while most of the lines above 518 nm are NH2 ([NI] 519.8 nm 
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Figure 2: The observed CN light curve in arbitrary units. In this figure and in the following 
ones, t is the time elapsed since the impact. The superimposed dashed lines are examples of 
model fitting. In all cases v = 400 m/s and = 0.2 d _1 . Model 1 (red) is computed with a v 
= 100 m/s, a p = 2 d" 1 and a g /a p = 4. Model 2 (green) is computed with a v = 100 m/s, a p 
= 0.1 d _1 and a g /a p = 100. Model 3 (blue) is computed with a v = 200 m/s, a p = 5 d _1 and 
a g /a p = 4. Model 4 (dotted pink) has no grandparent (a v = 200 m/s, a p = 0.1 d _1 ) and is 
clearly unable to reproduce the data. Fig. [2b shows the better agreement of model 1 with the 
two later UVES measurements (scaled to the HIRES slit area). 
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Figure 3: The observed OH light curve in arbitrary units. The triangles are scaled data from 
Schulz et al. (2006). Model 1 (red) is computed with v = 400 m/s, a v = 200 m/s, = 0.4 d _1 , a p 
= 0.5 d _1 and a g /a p = 40. Model OM1 (blue) uses the same parameters but it is adapted to the 
circular aperture of 6" (3900 km) used by Schulz et al. (2006) and scaled to correspond to model 
1 at large t (with a circular aperture, Eq. 5 must be replaced by n d = rid (1 — \A — (tw/t) 2 ))- 
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Figure 4: The observed NH2 light curve in arbitrary units. The illustrated model is computed 
with v = 500 m/s, a v = 100 m/s, ad = 0.6 d _1 , a p = 10 d _1 and a g /a p = 0.05. 
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Figure 5: The observed NH light curve in arbitrary units. The illustrated model is computed 
with v = 550 m/s, a v = 150 m/s, ad = 0.8 d _1 , a p = 2 d _1 and a g /a p = 2. 
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Figure 6: The observed C3 light curve in arbitrary units. The illustrated model is computed 
with v = 400 m/s, a v = 200 m/s, ad = 0.3 d _1 , a p = 10 d _1 and a g /a p = 0.01. 
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Figure 7: The observed C2 light curve in arbitrary units. The illustrated model is computed 
with v = 650 m/s, a v = 200 m/s, ad = 0.3 d _1 , a p = 1 d _1 and a g /a p = 4. 
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Figure 8: The observed CH light curve in arbitrary units. The illustrated model is computed 
with v = 400 m/s, a v = 200 m/s, ad = 8.0 d , a p = 2 d _1 and a g /a p = 15. 



n | i i i i | i i i r 



X 

2 2 



n i i i | r 
[01] 



1 - 



_l I I I I I I I I I I I I I I I I L_ 



0.00 0.05 0.10 0.15 

t (days) 




*.00 0.20 



0.40 0.60 0.80 

t (days) 



Figure 9: The observed [01] 5577 A light curve in arbitrary units. The illustrated model is 
computed with v = 400 m/s, a v = 100 m/s, ad = 0.4 d , a p = 4 d _1 and a g /a p = 10. 
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Figure 10: The observed dust light curve in arbitrary units. 
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Figure 11: Phase diagrams of the various species. The phase origin coincides with the impact, 
and the period is 1.709 d. The measurements obtained during the first post-impact nights (first 
night at phase ~ 0.45, second night at ~ 0.05) are denoted by circles. The curve is a fit with 
a three-source model optimized for the whole data set, not for a particular molecule (see text). 
Hence some fits may look rather poor, because the S/N of the data is bad (e.g., OH). The dust 
does not show periodic variations, but the June data appear to be systematically higher. The 
latter are denoted by a different symbol (x). The flux has been scaled to reach ~ 1 at phase 



origin. 



34 



